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Abstract 

It is shown that a set of functions which characterise the Lax hierarchy of non-linear 
equations may be represented in terms of the eigenstates of the potential which satisfies 
the generalised KdV equation. Such a representation leads to sumrules relating integrals 
involving the soliton potential and its various derivatives to sums involving the boundstate 
eigenvalues of the Schroedinger equation for the reflectionless potential. A new hierarchy of 
functions, which is in a sense dual to the Lax hierarchy, is identified. It is shown that time 
dependent equations involving the dual functions may be established which permit solutions 
related to an iV-soliton structure similar to that for the Lax hierarchy but with a different 
'speed' for the solitons. 



1 Introduction 

The Korteweg-deVries (KdV) equation was first discovered in the study of water waves. The 
KdV and the related equations with higher order nonlinearity, which are members of the Lax 
hierarchy (Lax 1968), have played a fundamental role in the study of nonlinear systems because 
they simulate many physical systems (Scott et al 1973), admit many conservation laws and the 
multi-soliton solutions can be given in analytic form. The KdV equation and its generalisation 
Kadomtsev-Petviashvli (KP) equation have also played an important role in pure mathematics 
because of their connection to algebraic curves, Jacobian varieties, vector bundles on curves, 
Schur polynomials and infinite dimensional Grassmannians (Mulase 1984). 

The connection between N-soliton solutions of the KdV equation and reflectionless potentials 
with TV boundstates in non-relativistic Quantum Mechanics is well known (Kay and Moses 1956, 
Gardner et al 1967, Scott et al 1973). Eventhough the reflectionless potentials belong to a 
rather restrictive category, they have nevertheless proved to be useful in some areas of Physics. 
The property of asymptotic freedom in Quantum Chromo Dynamics guarantees that quarks can 
not exist freely suggesting that inter-quark potentials are confining potentials and therfore have 
vanishing reflection coefficients. The zero angular momentum s-states of a confining potential 
in 3-dimension may be viewed as the odd states of a symmetric reflectionless potential. Such a 
picture has been used by Thacker et al (1978), Quigg et al (1980) and Quigg and Rossner (1981) 
to use experimental data for quarkonium systems as input information for the iV-soliton KdV 
algorithm to construct quark-antiquark potentials. 
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We consider higher order KdV equations given in the form (Lax 1968, Sawada and Kotera 1974, 
Caudrey et al 1976) 

dt m dx ^ 

where [Lj] satisfy 

t =(£-«'S-€) V.. » = (2) 

and i m is the 'time' parameter of the m th member of the hierarchy For example m = 1 leads 
to 

r) 2 JJ 

(3) 

and the third order KdV equation in the form 

dU d 3 U ^dU 

The solutions to the non-linear equations of the KdV and other members of the Lax hierarchy 
may be used as potentials in linear Schrocdingcr equations and their spectral properties may be 
studied by solving (using units in which mass = 1/2 and h = 1) 

H = +U , H 4, k = E k V> fe . (5) 

The N-soliton solution of eq. (1) may be viewed as a reflectionless potential U which supports 
N boundstates of the Hamiltonian operator H. Lax (1968) has shown that when the 'time' 
evolution of U is governed by eq. (1), if the 'time' evolution of the eigenstates tp k is governed 

by 

B m ipk = i-KT- i ipk{x,t) = exp(-iB m t) ip k {x,0) (6) 
where the operator B m satisfies the commutator relation 

then the eigenvalues E k of H are independent of t m and the eigenstates remain normalised, but 
the normalisation constants of the boundstates and the reflection and transmission coefficients 
for positive energies acquire a t m dependence. B m is a Hermitian operator which may be inter- 
preted as the generator of 'time' evolution which propagates the potential U according to the 
KdV equation or another higher order non-linear equation arising from eqs. (1) and (2). This 
propagation of U by B m is distinct from the usual time evolution of the Schroedinger eigenstates 
by the Hamiltonian H which propagates particles through a fixed potential. The hermiticity of 
B m ensures unitary 'time' evolution of the eigenstates. 

For the case m = 1, which leads to the third order KdV, the explicit form of B\ is given by 

*-'(-^^s + 'i)- < 8 » 

If the potential evolves in t\ according to eq. (4) and the eigenstates evolve in t\ according 
to eqs. (6) and (8) then the eigenvalues [Ek = —7^] in eq. (5) are independent of t\, but the 
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normalisation constants [Ck] which determine the behaviour of [ipk] as x — > ±00 and the reflection 
coefficient R for positive energies E — k 2 depend on t\ as given by (Scott et al 1973) 



CkilkM) = C k (j k ,0) exp(-4 7 |i 1 ), R(k,h) = R(k,0) cxp(8ik 3 t 1 ) . (9) 



Similar results hold for other members of the Lax hierarchy. One of the aims of this paper is to 
elucidate the explicit forms of [B m ] and [L m ] for all members of the Lax hierarchy. 

The plan of the paper is as follows : In section 2 of this paper we discuss some interesting 
properties of the well known N soliton solutions of the Lax hierarchy using a simple self-contained 
approach. In section 3 we find a representation of [L m ] of the Lax hierarchy in terms of the 
eigenstates of the N-soliton potential. We show that such a representation leads to sumrules 
involving the eigenvalues of the soliton potential. In section 4 we study the 'time' development 
of the soliton solutions and find an explicit expression for the eigenstate time evolution operator 
B m for the entire Lax hierarchy. The method used for deriving the time dependent equation 
suggests that it is possible to generalise the Lax hierarchy further. In section 5 we consider the 
concept of a dual to the Lax hierarchy and study the soliton structure of the first member of the 
new hierarchy. Section 6 contains a discussion of the main results of this paper. 

In this paper the word 'time' and the symbol t stand for a parameter which appears in the 
description of the evolution of the potential and does not refer to the coordinate t which is 
canonically conjugate to the energy in Quantum Mechanics. With this clarification we drop the 
quotation marks on 'time' in the rest of the paper. 

2 N-soliton solutions of the Lax hierarchy 

The N-soliton solutions of the KdV equation and all other members of the Lax hierarchy can 
be given in many equivalent forms. At any given value of a parameter t m , which is generally 
referred to as 'time', the 7V-soliton solution may be viewed as a reflectionless potential and the 
eigenstates of the Schroedinger equation of the potential may be studied. Here we consider a 
representation of the reflectionless potential with N boundstates used by Thacker et al (1978). 
The procedure for the construction of of reflectionless potentials with boundstates is a special 
case of a general procedure for finding a new potential by adding boundstates to a given potential 
(Sukumar 1986, 1987, Baye 1987) and the time evolution of the general potential will be studied 
in a future publication. In this paper we specialise to the case of the reflectionless potential. 
We suppress the dependence on t m for now and include it explicitly when studying the time 
evolution later in this paper. 

Eventhough the reflectionless potentials have been studied extensively in the literature we outline 
the key steps of the construction here to elucidate the connection to the generalisation of this 
construction in section 4. Starting from a potential Uq = with no boundstates, (viz) a free 
particle, for which the reflection coefficient vanishes for all positive energies, the solutions to the 
Schroedinger equation for energies [Ek] = [—7^], given by 



C k exp(--f k x), k = 1,2, N 



(10) 



may be used to to define a matrix whose elements are 




\ k {x)\i{x) 
Ik + li 



k,l — 1,2,.., 



N . 



(11) 
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>The solutions [ipi(x)] to the system of linear equations 



N 

Y^Au V; = A fc , fc = l,2,....,JV (12) 
1=1 

may be constructed and it can be shown that by using Kramer's rule to express ipi as the ratio of 
two determinants, (viz), the ratio of the determinant obtained by replacing the I th column of A 
by Afc to the determinant of A. Using the property that the differential of A^i is — A^A; and that 
the differential of a determinant is a sum over the N determinants constructed by differentiating 
each column in turn it can be shown that 

N N N N N „ . „ 

W ( X ) = -E^^ = -EE^[^- 1 L^ = EE^- 1 L*^ L = £ Ind ^- ( 13 ) 

1=1 1=1 k=l 1=1 k=l 

•It may be shown that [ipi] define the N boundstates of the Schroedinger equation for a potential 
U given by 

dW d 2 

U( X ) = -2— = -2- IndetA . (14) 
To prove this eq. (12) can be considered in the form 

N AfeA z 

i>k = A fe - E — -, — i't ( 15 ) 
t~l lk + 11 

and differentiated twice and rearranged to give 



Using eqs. (11)-(15) this can be brought to the form 

A *i(-^2 +li+U)^i(x) = (17) 
which can be used to infer that if det A ^ then ipi must satisfy 

showing that is a solution in the potential J7 for the energy — jf. 

• Eqs. (11) and (12) may be used to show that [tpi] satisfy the boundary conditions appropriate 
for boundstates at x — > ±oo and hence [tpi(x)] are the N boundstates of the potential U(x) defined 
by eq. ( 14 )■ In the limit 

x -> oo , A M -> 4 ( , dct A ->• 1, Vfe -» A fe = C fc cxp(-7 fc a;) -> (19) 

and it may be shown that in the limit 

^ ^ d — (nf ) (nn(^) 2 ) . <»> 



4 



A similar evaluation of the determinant in the numerator of the expression for ipk expressed as 
the ratio of of 2 determinants (solution to eq. (11) by the application of Kramer's rule) may be 
carried out to give 

~-.*(^(n^)G&)-°- <»> 

The choice of the normalisation constants Ck given in Thacker et al (1978) 

C l = TT (li+lk) (m 
2 lk 11 |(7/ ~lk)\ K ' 

corresponds to the requirement that \ipk{R)\ = |V'fe( — R)\ m the limit R — > oo, which in turn 
corresponds to the condition that U (x) be a symmetric potential. With this choice a symmetric 
reflcctionless potential with N boundstates may be constructed. It has been shown by Thacker 
et al that [ipi] defined by eq. (12) are indeed normalised to unity in the interval [— oo, oo]. Here 
we note that since the potential defined by eq. (14) depends only on the second derivative of the 
logarithm of a determinant, two determinants which defer by a factor exp(ax + b) will lead to the 
same potential whatever be the values of the constants a and b. Hence there are many equivalent 
representations of the matrix A which lead to the same potential U (Sukumar 1986). 

•It may be shown that the reflectionless potential with N boundstates may be represented as a 
weighted sum over the N boundstate probability densities. To establish this eq. (15) may be used 
to get 

N N N j IV IV 

^ifcV'fe^fe - Afe) = -j2J2^ kXk ^ x 3~r — = -sEE^^^ = — 2~ ' ^ 

fe=i fe=ij=i Ij-rlk fe=ij=i 

Differentiation of eq. (12), multiplication from the left by tpk and summation over k yields 



N N „ . N N N „. 

EE** Ak < -t - E** Xk E*< * = j> *± . (24) 

fe=l 1 = 1 ' k=l 1=1 k=l 



The first term can be simplified to a sum over A; multiplying the derivative of ipi and rearranged 
to give 

dipt . d\i 



i=i 

The definition of W in eq. (13) and U in eq. (14) can then be used to show that 

E( A <^ + *» = = 2 ■ (26) 



Thus we can establish that 



E A '^ = T + — ' ( 2? ) 

E^ = "E^'* = — — . (28) 
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Eqs. (23) and (28) may now be used to derive the expression 

N 8W 8 / N \ 8 2 

U = -4j># = -2— = 2-(]>>A fe ) = -2— IndetA (29) 

fe=i fe=i 
exhibiting the various forms in which the potential may be expressed. 

• The expression ofU as a weighted sum over boundstate probability densities leads to a sumrule. 
It is clear from the asymptotic limits given in eqs. (19) and (20) that 

£**->«, W = -Lt x ^ x C 2 k cxp(-2 lk x) -> (30) 

k 

8 N 
£«*->-«> W = Lt x ^_oo —In det A = -2^ 7fe . (31) 

X k=l 

Integration of eq. (29) then gives the sumrule 

f°° / \ N 

/ U(x) dx = -2 (W(oo) - W{-oo)) = -4 ^ 7fc (32) 

fe=i 

W /-OO 

= -4 E> / ^(*) • ( 33 ) 
fe=i J -° 

Eq. (32) shows that the area enclosed between the cc-axis and a graph of the potential U (x) is 
simply related to the binding energies of the boundstates of the potential U. Eqs. (32) and (33) 
taken together confirm that the eigenstatcs [ipk] defined by eq. (11) are normalised to unity in 
the interval [— oo, +oo]. 

Since U{x) is the lowest member of the [Lj] hierarchy, (i.e) L = U, it may be possible to find a 
representation of the higher members of the [Lj] hierarchy in terms of the boundstate probability 
densities. We examine this possibility next. 



3 Lax hierarchy and the eigenstates of U 

We now consider the differential equation satisfied by the probability density associated with the 
normalised boundstate ^fc at energy E k = —J k , for any potential U : 

P k (x) ee (Mx)) 2 , j^Vfc = (U + l 2 Mk ■ (34) 
It can be established that 

= 2UH + 2^2i\. + 2 



Qx2 - 2UP k + 2Y k Pk + 2^j . (35) 

By taking a further derivative of eq. (35) a differential equation satisfied by P k can be established 
in the form 

(£ - < - €) * - ^ t ■ <-» 



G 



It may be verified that the three linearly independent solutions of this third order differential 
equation are 4>%,ipl and V'fcV'fc where tpk is the second linearly independent solution to the dif- 
ferential equation satisfied by ipk- The weighted sums over the probability densities defined 

by 



N 



Qj = -4^/3 fc (2 7fc ) 3 P k 



where are arbitrary functions of k, therefore, satisfy 

d 3 .„ d „dU\ „ dQ j+1 



) Qj = 



dx 



(37) 



(38) 



which is identical to the differential equation linking the different members of the [L] hierarchy 
(eq. (2)). For the particular choice = [ 7 fc], we find using eq. (29) that 



Qo 



U 



(39) 



and the [Q] hierarchy defined by eqs. (37) and (38) becomes identical to the [L] hierarchy defined 
by eq. (2). 

• [Lj] may be represented in terms of the boundstate probability densities in the form 



N 



Lj{x) 



which leads to the sumrules 



k=i 



N 



/oo lv 
Lj(x)dx = -2^(2 7fc ) 2j+1 , j = 0,l,2,... 
-°° fe=i 



(40) 



(41) 



The first 3 members of the [L] hierarchy obtained by iteration of eq. (2) are explicitly given 

by 



U 
Li 

Li 



U 

d 2 U 
dx 2 
d 4 U 
dx 4 



-3U 2 



-0 - + - 



lead to the sumrules 



J — ( 

/OO 



N 



U dx 



U 2 dx 



dU\2 



dx 



+ 10C/-' 



-4 £ 7 fc 
k=i 

16 fit 

k=i 

N 

dx = -64 ]T 7fe 5 



(42) 
(43) 

(44) 

(45) 
(46) 
(47) 



fc=i 
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It is easy to check that these sumrules hold good for the 1-soliton potential U = — 2^1sec\? r )\x 
for which all the integrals can be carried out analytically. 



The procedure given above may be used to establish N sumrules from the first N members of 
the [L] hierarchy and if the reflectionless potential with TV boundstates is known then the first N 
sumrules may be used to determine the N binding energies from the various integrals involving 
U and its derivatives. This is yet another striking property of the soliton solutions of the Lax 
hierarchy. 

4 Time evolution of the potential and eigenstates 
4.1 Evolution of the potential 

The potential U and the boundstate solutions considered in sections 2.1 and 2.2 are time inde- 
pendent because the basis functions in eq. (10) have no time dependence. If we consider basis 
functions X k which depend on a parameter t in the manner given by 



then eqs. (11)-(18) acquire a dependence on t. All and A appearing in sections 2 and 3 may be 
viewed as corresponding to functions evaluated at t = and the time parameter can be restored 
by replacing all functions of ip(x) and A (or) by ip(x,t) and X(x,t). Eqs. (11) and (48) may be 
used to find the i-derivative of A k i in the form 



The i-evolution of ij)k{x, t) may be studied by differentiating eq. (12) with respect to t and using 
eq. (49). These algebraic manipulations lead to 




(48) 



d_ 

di 



Am = (a k + at) (A k i - S k i) 



(49) 




(50) 



The time evolution of ipi can now be given in the form 




(51) 



fe=i 



Multiplication of cq. (50) by ip k and summation over k yields 



N N 



din 
dt 



N 



+ ai ipi J = a k V>fc 



k=l 1=1 



fe=l 



which can be simplified using eqs. (12) and (48) to give 



Using eq. (14) we can then establish the evolution equation for U in the form 

^= 4 ^E^ (53) 
fe=i 

• We note that this evolution equation is valid for arbitrary choices of [a k ] and in this method of 
construction of U the boundstate eigenvalues of U are independent of t which implies that this 
time evolution ofU can exist alongside a unitary evolution of the eigenstates ofU. The potential 
constructed using eqs. (11)- (14) with the basis functions given by eq. (4-8) can now be viewed as 
an N-soliton potential. The soliton structure ofU(x, t) can be extracted using the procedure used 
by Gardner et al (1914)- The non-linearity of the evolution ofU is implicit in eq. (53). However 
for special choices of [a k ] the non-linear character can be made explicit. The consequences of 
this distinction will be further elaborated in section 5 of this paper. 

For the special choice of [a k = 2 2m 1 2m+1 ], m = 1, 2, .., the sum in the right hand side of eq. (53) 
can be identified as L m , a member of the Lax hierarchy defined by eq. (40). We now use the 
symbol t m to identify this particular t-evolution. Comparison with eqs. (40) and (53) then shows 
that the time evolution equation for U(x,t m ) is now given by eq. (1). For all members of the 
Lax hierarchy eq. (2) may be used to express L m in terms of U and its spatial derivatives and 
the non-linear character of the evolution equation (1) can be made explicit. The time dependent 
[L m ] can be given in the form 

T , . , , 2 . . „ dW „ d ^ , . . , . . „ d d , , . 

L m (x,t m ) = -A2_^a k ip k {x,t m ) = 2— = -2— ^\ k (x,t m )ip k (x,t m ) = 2 — — IndetA 
_^ ot m dt m dt m ox 

(54) 

which are all equivalent and provide the generalisation of eq. (29). For the case m = the 
evolution given by eq. (1) is simple since in this case a k — j k and L — U and U and W are 
functions of {x — i). 

For integer values of m, the potential U(x,t m ) is a solution of the m th member of the Lax 
hierarchy defined by eqs. (1) and (2), which is also referred to as the (2m + l) th order KdV 
equation. We note that the limiting values of det A in the asymptotic region x ±oo given by 
eqs. (19) and (20) are still valid when the t m dependence of A is included. Hence using eqs. (19) 
and (20) for the asymptotic limits of det A and eq. (48) to evaluate the time derivative of In Afe 
it can be shown that 

N 

L m (x,t) dx = -2^(2 7fc ) 2m+1 (55) 
fe=i 

which is in agreement with eq. (41), the sumrule derived in section 3. Since the eigenstates [tp k ] 
in eq. (40) are normalised to unity for all t m , the sumrules given by eq. (41) are independent of 
t m and are valid for any value of t m . 



J — : 



4.2 Evolution of the eigenstates 

We next examine the time evolution of the eigenstates of the time dependent potential U (x, t m ) 
when the basis states Afc evolve according to eq. (48). The explicit form of the time evolution 
for the third order KdV is given by eqs. (6) and (8) which can be simplified using eq. (18) to 
the form 



M^h) = -4 7fc — + [2L — - —4> k ) . (56) 
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We now show how equation (56) arises from eq. (51) for the case m = 1 for which [a k = 4 7 j?] 
and find the form of the equation analogous to eq. (56) for other values of m. 

The procedure for adding TV boundstates to a potential Uq = may be reversed so that starting 
from a potential with iV boundstates in U we can find a potential with no boundstates in the 
form 

d 2 

U = U - 2 — IndctB = (57) 



B k i = s k i + / tpk(y) My) d y = / i>k{y) i>i{v) d v ( 58 ) 

and the two sets of functions A and tp are now related by 

N N 

1=1 fc=i 

so that 

JV N N N N „ „ 

1=1 1=1 fe=i /=i fc=i 

Comparison of the two sets of equations corresponding to the addition and the removal bound- 
states shows that the matrix [B] is the inverse of the matrix [A]. Hence 

[A-\ k - B lk = f My)My)dy (6i) 

J — 00 

which can be further simplified using the Wronskian relation 

^ it ~ 7& ^ fe = (7 ' " ^ ^ d » • fc ^ • (62) 

The steps taken for establishing the time evolution equation (51) starting from eq. (48) may 
be repeated to find the x derivative of tp, by the replacement a k — > — 7 fc in eq. (51), in the 
form 

^) = - 2 E[^ _1 ]/fc 1* ^ ( 63 ) 



{dipt 

fe=i 



For the KdV-Lax hierarchy [a^] may be given in the form 

a k = A m 7 2m+1 = 4 m lk ( 7 2m - 7 2m ) + 4 m 7fc 7 2m (64) 

and a factor ( 7 | — 7 2 ) may be extracted from the first term on the right hand side of eq. (64) so 
that the combination [^4 _1 ] ;A , ot k in eq. (51) may be expressed in terms of a Wronskian. These 
algebraic manipulations lead to 



dt 



N 

1+^^=4™ 7i 2m (2E [ A ^]ik ^ A) +A t (65) 



fe=i 

N 



a ; = 2^ e d -t ~ d t **) (E^ -2 ^ m ' 23 h ^ • (66) 
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Using eq. (63) for the x-derivative of V and eq. (40) for [L k ] it is possible to simplify the 
evolution equation for ip to the form 



a*, _ _ {2iir g + E (2 7i)2 ,-, ( - ^ „, + 2 ^ f ) . (67 , 



at 



It is readily verified that for m = 1 this expression reduces to eq. (56) for the third order 
KdV. 

• The time evolution of the eigenstates ofU(x,t m ) for the entire KdV hierarchy is given by 

Q 

-^—ipk(x,t m ) = -i B rn fa(x,t m ) 
ut m 

B m fa(x,t m ) = * ( 2 7*) 2m ^ + * E ( 2 ^) 2j_2 ( 2L ™-^ ^g^) ■ (68) 

We have established this equation from first principles. It is possible to verify directly that this 
is consistent with eqs. (5) and (7) as follows. We can use eq. (5) to show that 

[B m ,H] fa = + 7fc) B m fa . (69) 

The right hand side of eq. (68) consists of three types of terms and the effect of (H + 7^) acting 
on each of the terms can be found and simplified using eq. (5). Eq. (2) for the [L] hierarchy 
can then be used to futher simplify the expressions. These algebraic manipulations lead to the 
result that 

BT 

[B m ,H] fa = -i-^r V>* • ( 7 °) 

We can also use eqs. (29), (40) and (68) to show that 
dU A d , 2 



+ 2 



dx ^ \ dx J J dx 

It is easy to show that all the terms in the sum cancel each other when the sum is expanded. 
Thus we can verify that 

dT 811 
[B m ,H] fa = -i-^fa = i — fa (72) 

showing that the time evolution of fa given by eq. (68) is consistent with eqs. (1) and (5) 
and (7). The time evolution operator for the eigenstates of the iV-soliton potential of the Lax 
hierarchy can also be given in the form 

«.--«£(. 4*r + . t - %*) < - <*)'-' ■ m 
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5 A dual Lax hierchy for reflect ionless potentials 



We now examine the general evolution equations (48) and (53) again. It was noted earlier that the 
m th member of the KdV hierarchy arises from the choice of \a.k = 2( 2 " l )7^ 2 " l+1 - ) ] . As emphasised 
earlier, eq. (53) is valid for arbitrary values of [a^]. This observation allows us to consider an 
alternate choice of [a^ = 2^ 2m ^^ 2m ^] leading to a hierarchy of functions L m ,m = 1,2,..., 
which correspond to the negative values of the index m of [L m ] considered earlier. Thus the 
members of the hierarchy defined by 



satisfy 



^ = (ifs - 2^)Z m , m=l,2,.. (75) 

ox V ox 6 ox Ox I 



and lead to the time evolution equation for the potential of the form 

dU dl m 



dt dx 



(76) 



• This hierarchy may be viewed as a dual hierarchy to the usual Lax hierarchy. In the Lax hierarchy 
starting from Lq — U higher members are found using eq. (2) by differentiation and integration. 
In the hierarchy defined by eq. (75) starting from Lq = U higher members are found by solving 
differential equations. In this respect the processes for generating higher members of the sequence 
for the two hierarchies differ in a fundamental way. The potential U is still represented by a 
weighted sum over the probability density of the normalised eigenstates. The method used by 
Gardner et al (1974) to exhibit the soliton structure of U(x, t) in the context of KdV still applies 
and it can be shown that in the asymptotic domain, t — > ±oo, x — > ±oo, the potential becomes 
a superposition of N localised solitons, localised around different regions of x, just as for the 
soliton solutions of the KdV hierarchy. 

We now consider the m = 1 member of this family, which can be shown to have a special 
significance since the non-linearity of the evolution equation for this case may be unravelled. 
The relevant equations for this case are: 

X k {x,t) = C k exp(-7 fc ir + -^-) (77) 



N 



, , ,x x A fc (x,t) Xi(x,t) 

A k i(x,t) = bki + ■ (79) 

Ik + li 



N 



^A k ii>i{x,t) = X k {x,t) , fc = l,2,..,JV (80) 
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from which it is possible to construct 

U(x,t) = -2 —\ndctA(x,t) = -4^7j ipj(x,t) (81) 

du dU 

~dt ~dx~ (83) 

The eigenvalues [— 7? ] of the Schroedinger equation for the potential {/ (x, t) defined by the above 
equations do not depend on time showing that the time evolution of [tpj(x,t)] implicit in eqs. 
(81)-(83) is a unitary evolution. 

The evolution equation for the eigenstates can be studied by starting from eq. (51) and following 
steps similar to the steps taken for the Lax-KdV hierarchy (eqs. (61) - (66)). It can be shown 
that 

dipt 1 dipi 1 /dLx dipi\ 

~dt = ~Wf -dx~ + 4^fU^ - 2Ll ^) (84) 

which is the analogue of eq. (56) for the eigenstate evolution when U satisfies the KdV equa- 
tion. 

The non-linearity implicit in the evolution equation for U given by eq. (83) can be made more 
explicit by the following reasoning. Since L\ also satisfies 

(dx 3 ^ dx dx) _1 dx dx 

it is evident that 

We now observe that the solutions to the Schroedinger equation for the potential U(x, t) at zero 
energy given by 

= U{x,mx,t) (87) 

may be used to show that 

Comparison of eqs. (86) and (88) shows that a solution £ to eq. (87) which satisfies suitable 
boundary conditions may be found such that 

u = «^ . m 

Eqs. (83) and (87) can now be combined to show that 

dt\i dx 2 ) dt dx ^ dx' v ' 
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Thus we have shown that the nonlinear equation 



9a; 2 dt dt dx 2 dx 

has solutions of the form 

i 2 M = 1 2 £ ^M. (92) 

The potential C/ is a symmetric function of x at i = if Ck are chosen according to eq. (78). 
The structure of U, [ipk] and £ in the aymptotic domain, t — > ±oo, a; -4- ±oo, may be analysed 
by using the same method as Gardner et al (1974) and Thacker et al (1978). In the asymptotic 
domain the potential U{x,t) is a superposition of N shifted soliton solutions of the form 

u k ->• -27 2 sech 2 (7 fc y fc =F 4) (93) 
t 



Vk 



= (x - — j) , fc = l,2,..,JV (94) 

< k 



k-1 , , N 



localised around the N different regions of x where a particular y^ is finite while all the other 
[yj — > ±oo, j 7^ fc]. The normalised eigenstates [tpk] are also localised around the same distinct 
regions of x and have the form 

V>fc -> y^sech( 7fe2/fe t 4), fc = l,2,..,iV. (96) 

Hence 

£ 2 -> tanh 2 ( 7fe y fe =F 4) , fc = l,2,..,iV (97) 

(i.e) £ 2 is also localised around N different regions of x where [yu] are finite. Thus the non-linear 
equation we have derived has solutions with a structure which is related to a soliton structure 
similar to that for the KdV hierarchy of equations. However for t — > ±oo the solitons have speeds 
different from those of the KdV solitons because of the altered time dependence carried from eq. 
(77) to eq. (94). 

£(x, t) is also the solution to the Schroedinger equation at energy E = for the potential given by 
eq. (81) which satisfies appropriate boundary conditions which are compatible with those satisfied 
by [i[)j(x,i)], as required by eq. (92). The zero energy solution may be directly calculated as 
follows. For positive energies E = k 2 it may be shown that eq. (15) may be extended to the 
form 

<P(E,x) = exp( ±J fcx) - J2 XliX) e ^ (± ' kx) M*) (98) 

and the steps similar to the passage from eq. (15) to eq. (17) may be performed to show that 
the tp(E,x) so defined is indeed a solution of the Schroedinger equation for energy E. We can 
now find the zero energy solution by setting k = in eq. (98) to get 

= m*) = i - E MM) ^ (M) - (99) 
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By squaring this expression and using eq. (15) it may verified that eq. (99) is compatible with 
eq. (92). £ satisfies the boundary condition Lt^^ £(x) — > 1 and using eq. (21) it may be shown 
that Ltx^-oo £(x) — > (— ) N . These boundary conditions are compatible with eq. (92). Thus 
£ may be expressed in cither of the two forms given in eqs. (92) and (99). The non-linearity 
of the time evolution of £, as shown by eq. (91), may thus be used to characterise the implicit 
non-linearity in the time evolution of U(x,t) through eq. (90). 

For the 1-soliton case the solution to the non-linear equation (91) is 

£ = tanh (71 yi ) , yi = (x - ^) ( 10 °) 

and using 

01 = V^Ti ex P ( - 7i Vi) j "01 = sech(7i j/i) , Ui = -2jf scch 2 (7i 2/1) (101) 

eqs. (86) - (92) and eq. (99) may be verified. 

For the 2-soliton case the solution to eq. (91) may be shown to be 

^ _ 72 sinh (71 yi) sinh (72 y 2 ) - 71 cosh (71 2/1) cosh (72 y 2 ) 
72 cosh (71 yi) cosh (72 2/2) - 7i sinh (71 j/i) sinh (72 y 2 ) 

yj = (* - 4^2) , i = 1,2. (103) 

and the corresponding potential and eigenstates for the case 72 > 71 are 



U 2 = - 2 W2 n2 7i; hi cosh 2 (71 yi) + 7 2 sinh 2 ( 72 y 2 ) ) (104) 



D 2 

D = (j 2 cosh (712/1) cosh (722/2) - 7i sinh (71 2/1) sinh (722/2)) (105) 

/ FF /T~2 2T sinh (72 2/2) 

V-i = -Jy V(^2 "7i 2 ) ^ (106) 



^2 



Vf - (107) 



For general TV, £ may be constructed using eqs. (99), (80) and (77) and may be shown to be 
related to a Schur polynomial of order N. 

So far we have examined the stucture of the non-linearity for the m = 1 member of the system of 
eqs. (74) - (76) in detail. It is clear that other members of this hierarchy also admit an iV-soliton 
structure for U and a unitary evolution of the eigenstates of U. The time evolution equation for 
the eigenstates of all members of the dual hierarchy can be shown to be 

f - - £ - £ ,r-' ( - * + > £) . doe, 

which is the generalisation of eq. (67) for the dual hierarchy. The dual hierarchy [L m ] discussed 
in this paper may be viewed as complementary to the Lax hierarchy [L m ]. 



15 



6 Discussion 



In this paper we have shown that for the KdV-Lax hierarchy of non-linear equations governing 
the time evolution of the potential U, a representation of [L m ] in terms of the eigenstates of 
the potential may be found which enables the identification of the soliton structure of [L m ] for 
t — > ±00. Such a representation leads to certain sumrules involving all the boundstates of the 
potential U. We have shown that the integral of [L m ] over the entire spatial domain [—00, +00] 
is proportional to the sum over the energy of the boundstate raised to the power of m+|. With 
sufficient number of such sumrules it is possible to solve for the eigenvalues of the boundstates 
entirely in terms of the integrals of [L m ]. We have also given explicit expressions for [B m ], the 
time evolution operator of the eigenstates, in terms of [L m ] and the derivative operator for the 
entire Lax hierarchy. 

It is evident from eqs. (48), (10)-(14) and (29) that, independent of the specific choice of [ah], at 
t = there is a potential U(x,0) which supports N boundstates of the Schroedinger equation. 
For t — > ±00 equations similar to eqs. (93)-(95) may be used to show that the V-boundstatc 
potential separates into N separated packets of potentials, each with its own depth proportional 
to the binding energy, each supporting a single boundstate or soliton for all choices of [oik]- For 
different choices of [a^] the dependence of the speed of the soliton on the energy is different. For 
the KdV-Lax hierarchy [oik] is proportional to [7^ m+1 ] giving rise to solitons which move in the 
time region t — > ±00 with speeds proportional to a positive integer power m of the energies of 
the boundstates that the solitons represent. The non-linearity of the evolution equation for U is 
of order (2m + 1) for this case. 

We have identified a new hierarchy [L m ] which leads to a new type of non-linear time evolution 
equations for reflectionless potentials which retains the soliton structure of the Lax hierarchy 
and produces a unitary time evolution of the eigenstates of the Schroedinger equation for the 
potential. An explicit expression for the time evolution operator for the eigenstates of the 
potential may be found for this case also in terms of [L m ] and [L m ] may be expressed in terms of 
the eigenstates. For the dual Lax hierarchy identified in this paper [oik] is proportional to ["fl~ 2m ] 
giving rise to solitons with speeds proportional to a negative integer power —to of the energies. 
The evolution equation for U is now implicitly non-linear and for to = 1 the character of the 
non-linearity can be unravelled in the manner discussed in section 5. For this dual hierarchy the 
reflectionless potentials at time t = given by U(x, 0) may be used to determine £(x, 0) by solving 
eq. (87) at t — subject to the boundary condition that £ 2 — > 1 as x — > ±00 and this solution 
may be used to find £(x, t) by solving the non-linear evolution equation (91) and this £(x, t) may 
be used back again in eq. (87) to find U(x,t) and then [ipk(x,t)]. In the KdV-Lax hierarchy the 
potential U plays the fundamental role in the sense that all other [L m ] are expressed in terms 
of U and its derivatives through eq. (2). For the dual hierarchy the zero energy solution in the 
potential denoted by £, satisfying the boundary conditions that |£| — > 1 as x — > ±00, plays the 
fundamental role in the sense that £ determines all the members [L m ] as is evident from eqs. 
(75) and (89). It is for £, in the m = 1 case, that it has been possible to establish an evolution 
equation which is explicitly non-linear as a function of time and this is perhaps the case for all 
the other members of the dual hierarchy. The other members of the dual hierarchy with m > 1 
deserve further study . 

There are other choices of the time evolution functions [at] which lead to other interesting 
evolution equations for U and [ipk], which retain the soliton structure of the KdV with the 
speeds of the solitons having other kinds of dependence on the energy but still providing unitary 
evolution of [ipk]- These aspects will be the subject of a future study. 
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